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' Abstract. It is shown that the quasi-one-dimensional Bose-Einstein condensate 

, is experimentally accessible and rich in intriguing phenomena. We demonstrate 

' numerically and analytically the existence, stability, and perturbation-induced 

dynamics of all types of stationary states of the quasi-one-dimensional nonlinear 
I Schrodinger equation for both repulsive and attractive cases. Among our results are: 

the connection between stationary states and solitons; creation of vortices from such 
states; manipulation of such states with simple phase profiles; demonstration of the 
fragility of the condensate phase in response to shock; and a robust stabilization of the 

■ attractive Bose-Einstein condensate. 

o 
o 
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' 1. Introduction 



A family of Bose-Einstein condensates (BECs) has been created in a series of ground- 
breaking experiments [|I], |, ^ over the last five years. These condensates, made 
of trapped, dilute gases of ^^Na, ^^Rb, ^H, or ^Li, are well-modeled by the nonlinear 
Schrodinger equation (NLSE) [||, also known as the Gross-Pitaevskii equation 0, |^. 
However, they have yet to be experimentally investigated in quasi-one-dimension (quasi- 
ID.) We define the BEG to be in the quasi-lD regime when its transverse dimensions 
are on the order of its healing length and its longitudinal dimension is much longer than 
its transverse ones. In this case the ID limit of the 3D NLSE is appropriate, rather than 
true ID or 2D mean field theories IM|. 



The ID NLSE is expected to be the correct description of the quasi-lD BEG [p!0| , [11 
It is a ubiquitous equation. Among other natural phenomena, it models Bose-condensed 
photons |1^, disordered media |13], helical excitations of a vortex line |14|, and light 
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pulses in optical fibers In this latter context soliton solutions to the NLSE have 



been used to greatly increase the efficiency of communications . 

Because the BEC is a coherent matter- wave it has a natural analogy with light. 
This insight has led to light-inspired experiments, as for example four-wave mixing 
and both continuous-wave and pulsed atom lasers [0, 0. The creation of a quasi-lD 
BEC would allow the exploration of a great body of experimental and theoretical results 
from fiber optics [|16| . 

In addition to the analogy with fiber optics there are several other reasons to 
explore the quasi-lD BEC. Cold atom waveguides have recently been experimentally 
realized [^0|, This work suggests the manipulation of matter waves in quasi-lD and 
perhaps the eventual creation of matter- wave circuits, just as fiber optics have come to 
replace wires in certain applications. The quasi-lD regime of the BEC is qualitatively 
different from its 2D and 3D counterparts: in the attractive case collapse is prevented, 
while in both the repulsive and attractive cases solitons are stablized. 

In a pair of recent articles all stationary states of the ID NLSE under box and 



periodic boundary conditions ||T0|, |Tl| were analytically elucidated. This provides 
a starting point for the analysis of the quasi- ID BEC. In the present investigation 
we impose box or periodic boundary conditions in the longitudinal dimension and 
box boundary conditions in the transverse dimensions. These boundary conditions 
model the potential quasi- ID regime of a number of present experiments: the atom 
waveguide the prolate harmonic trap in which is formed a cigar-shaped 

BEC |23|, the toroidal trap and finally an oblate harmonic trap with a barrier 
formed in the middle either by a second spin state of the same atom PBl 12^ or by a 



laser ||2^, in which is formed a pancake-shaped BEC with the center removed. Periodic 
boundary conditions are the best ffist model for toroidal geometries; box boundary 
conditions are a good starting model for cigar-shaped geometries. 

Our approach is to find the stationary states of the NLSE under simple, analytically 
solvable boundary conditions. These give physical insight into the corresponding results 
for potentials which require numerical solution. Here we demonstrate the merits of this 
approach by the extension of ID into quasi-lD, 2D, and 3D. Elsewhere we have applied 
such insight to the harmonic potential pSf, among others. 

The paper is organized as follows. In section Q we review the analytic stationary 
solutions in ID. In section ^ we numerically illustrate these solutions and their properties 
and predict how they should extend into quasi-lD and higher dimensions. In section ^ 
we numerically verify these predictions and thereby demonstrate that the quasi- ID 
regime is experimentally accessible for both repulsive and attractive condensates. In 
section ^ we show how quasi- ID-type stationary states extend into 2D and 3D and the 
rich phenomena that result. Finally in section B| we give our conclusions. 
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2. Solitons and stationary states 



We consider two cases of tlie BEC: repulsive and attractive atomic pair interactions. 
For our purposes they differ by the sign of the nonhnear term in the NLSE. The former 
is the one that has received the most experimental attention. The latter is generally 
unstable in 3D but stable in ID ||30| . 

The NLSE in ID has a number of special properties which are described in 
the mathematical literature. It is integrable may be solved exactly by the 



Inverse Scattering Transform ||33| , Q and has a countably infinite number of conserved 
quantities [pSj |36[| . We have taken advantage of its special properties in finding its 
stationary solutions. We first describe all stationary states of the ID NLSE subject 
to box and periodic boundary conditions. We then make clear their connection with 
solitons and provide physical interpretations of their form. 



2.1. Mathematical form 

We begin with the NLSE that describes a BEC of atoms of mass M, confined in an 
external potential V{r): 

f~ 2M^' + 9 \^ir, t) I' +V{r} ] ^(f, t) = iWt ^(f, t) (1) 

where \ip{r,t)\'^ is the single particle density such that p{r,t) = N [^/^(r, the coupling 
constant g = Anfi^aN/M, and a is the scattering length for binary collisions between 
atoms. The case of repulsive interactions corresponds to a > 0; that of attractive 
interactions to a < 0. 

The characteristic length scale of variations in the condensate wavefunction is the 
healing length ^: 

^^i8n-p\a\)-'/' (2) 

where p is the mean particle density. The BEC is in the quasi- ID regime when the 
lengths of the tranverse y and z dimensions of the trap, Ly and L^, satisfy the following 
criteria: Ly,Lz ~ ^e// and Ly,Lz « L^. The former ensures that the condensate 
remains in the ground state in the two transverse dimensions while the latter ensures that 
longitudinal excitations are much lower in energy than possible transverse excitations. 
Under these conditions one may make an adiabatic separation of longitudinal and 
transverse variables ||10|, Equation (0) then reduces to: 

[-iM'dl± I fix) p +Vix) ] fix) = tdifix) (3) 

where all terms have been made dimensionless: fix) is the longitudinal portion of the 
wavefunction; t = ih/2M^'^)t we will call the natural time; the ± refer to the repulsive 
and attractive cases, respectively; Vix) = (2M^^//i^)y(x); and ^^ff is an effective 
healing length, discussed in section ^A\ . If we furthermore assume a stationary state: 

[-i^eff)'dl± I fix) P +Vix) ] fix) = p fix) (4) 
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where fx = {2M^^/h^)fi is a dimensionless chemical potential which is now an eigenvalue. 

All stationary solutions to equation (|^) may be written in terms of Jacobian elliptic 
functions |38|. The properties of such functions are reviewed elsewhere [10, 38, 39[ . 



There are five normalizable symmetry-breaking solution-types to the stationary NLSE 
with a constant potential on a finite interval. They are pictured in figure |I[ The 
stationary solutions appropriate to longitudinal box boundary conditions are: 

„, , ^ .2jK(m)x e- I \ 

f{x) = A sn( ^ ^ +6\m) (5) 
f, X ^ ,2jK(m)x r , ^ 

f{x) = A cn( ^ ^ +5\m) (6) 

J^x 

for the repulsive and attractive cases, respectively. Equation is shown in figure |I](a) 
and equation (|^) is shown as the anti-symmetric solution in figure |T](c). A is the 
amplitude, j — 1 is the number of nodes, K{m) is the complete Jacobian elliptic integral 
which is the quarter period of these functions, 5 is the offset, here equal to zero, and 
< m < 1 is the Jacobian elliptic parameter. The general notation sn(M | m) is 
standard for Jacobian elliptic functions [^, 

Under periodic boundary conditions equations @ and (|) are still solutions, but 
j is required to be even and the number of nodes is j, i.e. j = 2,4,6, ... Translational 
symmetry is restored, so the offset 6 is no longer required to be zero. This leads to a 
Kosterlitz-Thouless-type entropy in ID [T0[| . 

In addition to the above two solution-types there are three nodeless, symmetry- 
breaking solution-types under periodic boundary conditions. The first, a solution to the 
attractive case, is real: 

„, , ^ , ,2jK(m)x r I N /w\ 
/(x) = A dn( ^ ^ +S\m) (7) 

^x 

An example is shown as the symmetric solution in figure |l|(c). The other two types are 
intrinsically complex. f{x) = r{x) exp(z0(a;)) and for the repulsive and attractive cases, 
respectively: 

r{xf = A\l - 7 dn\^^^^^ + 6 \ m)) (8) 

Lx 

r{xY = A\dn\'^^^^^^ \m)-j{l-m)) (9) 
Lx 

where in both cases the phase must be found by numerical integration from the equation: 
= (10) 

The phase and amplitude are shown in figure |I](b) and [l|(d). In the repulsive case 
y4^7 is the depth of the density minima below the constant background. When 7 = 1 
equation is recovered. In the attractive case 7 interpolates between the real anti- 
symmetric and symmetric solutions in equations (^) and (|^. In both cases < 7 < 1. a 
is a constant of integration, j is the number of density minima or maxima, respectively. 

All five solution-types are found by solving equation (||) subject to normalization 
and boundary conditions, and are described in detail in references and |Tl[ . 
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2.2. Physical interpretation 

The nonlinear mean field term in equation (^) causes a spreading out or a clumping 
together of the wavefunction between nodes for the repulsive and attractive cases, 
respectively. Thus the stationary solutions to the NLSE with box boundary conditions 
are in one-to-one correspondence with the sinusoidal solutions to the particle-in-a-box 
problem in linear quantum mechanics. The number of nodes is j — 1 and the complete 
set of stationary states is described by equations © and where j G {1, 2, 3, ...}. 

Likewise, if j is restricted such that j G {2, 4, 6, ...}, equations (|^) and (H) are in one- 
to-one correspondence with the sinusoidal solutions to the particle-on-a-ring problem in 
linear quantum mechanics. As the ring has rotational symmetry these solutions are 
symmetry-breaking. We note that, as in the linear Schrodinger equation, there are also 
complex, constant-amplitude, plane-wave solutions on the ring. 

The Jacobian elliptic parameter m governs the strength of the nonlinearity. As 
m — > 0"'", sn ^ sin and cn cos, respectively. This is the linear, sinusoidal limit. 
As m ^ 1^, sn — > tanh and cn — > sech. These are the dark and bright soliton, 
stationary solutions to the NLSE on the infinite line, respectively [^S], This shows the 



connection between these stationary solutions and solitons. In the box the j*^ stationary 
state is a j — 1 or j soliton-train in the repulsive and attractive cases, while on the ring 
the 2^'*^^ stationary state is a 2j soliton-train. As we have shown elsewhere soliton- 
trains interact in a manner similiar to single solitons. And as with single solitons their 
form and phase relationships are quite robust. 

The other three solution-types have no analogue in the particle-on-a-ring problem in 
linear quantum mechanics. They exist only on the ring and are nodeless and symmetry- 
breaking. The complex solutions have a monotonically increasing phase. Each has a 
complex-conjugate, degenerate partner. The complex solutions for the attractive case, 
in particular, are an entirely new stationary solution-type for the NLSE. 

The complex solutions for the repulsive case are interpreted as density-notch solitons 
moving with speed c on the ring with an opposing momentum boost of the condensate 
of speed — c, which results in a stationary state in the lab frame. Density-notch solitons 
have a speed between zero and the Boguliubov sound speed, ranging from maximal to 
zero depth, respectively |4^. Those not of maximal depth are called grey solitons, while 



those which are of maximal depth and therefore form a node are called dark. Figure |l|(b) 
shows the bounded, quantized version of a grey 2-soliton-train. 

At typical experimental trap sizes the single density-notch stationary state on 
the ring is the lowest energy, symmetry-breaking excitation above the real, constant- 
amplitude ground state. When any of these stationary excited states is perturbed it 
gives rise to soliton- type motions |^ . Recent reports suggest that such motions can be 



induced in repulsive BECs by optical phase engineering of the condensate phase |2^, |4TI] 



All attractive symmetry-breaking, longitudinally periodic, stationary solution- 
types, i.e. the anti-symmetric and symmetric ones shown in figure [^(b) and the complex 
one shown in figure |l](d), are described by the Cj point symmetry group, where j is the 
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number of peaks. There are j nearly degenerate solutions. For even j there is a real 
symmetric-anti-symmetric pair and (j — 2)/2 degenerate complex pairs. For odd j there 
is a real symmetric solution and {j — l)/2 degenerate complex pairs. Thus by group 
theory these three solution-types form the complete set of stationary states made of 
evenly spaced peaks. 

As the attractive BEC has been made in 3D but not in quasi- ID, these stationary 
solutions are open to experimental investigation. 



2. 3. Phase engineering of dynamics 

There are two ways to impart velocity to a density-notch soliton. One may boost the 
constant background via a phase ramp and the notch will drift with the resulting super- 
current. One may also apply a phase jump and the resulting soliton speed with respect 
to the background will be c = Cmax sin(5/2), where Cmax is the Boguliubov sound speed 
and 6 is related to the phase difference across the moving density notch . Both types 



of phase profile are apparent in figure ^(b). The phase ramp may be seen in the constant 
background slope in the phase. The phase jump is shown over each of the two notches. 
The velocities cancel and a stationary state results. 

To impart a collective velocity to trains of density-notches one applies an equal 
phase difference across each of the member notches or a phase ramp across the whole 
train. One may also treat the members of the train individually by applying varying 
jumps and thereby study solitons interactions. A bright soliton may be set in motion by 
a phase ramp but not a phase jump. For a phase ramp exp{ikx) the resulting velocity 
c (X k. There is no upper bound on the velocity of bright solitons. As with density- 
notches, a train of bright solitons may be manipulated collectively by applying a phase 
ramp across the whole set or individually by applying different phase ramps across the 
members of the train. 

Soliton-trains are thus collective excitations which can be directly manipulated by 
simple phase profiles. One may manipulate the members of the train individually or 
collectively. This makes clear the connection between the above stationary states and 
solitons. 



3. One dimension 



3.1. Stability 

We first demonstrate that the five solutions types described in section |^ are indeed 
stationary. In figure |^ we propagate them in time numerically. In the numerical work 
presented for one, two, and three dimensional solutions of the NLSE the condensate 
wave function was expanded as ip{r,t) = T,iCi(t)(f)i{r). The 0j(r) were chosen as a 
co-ordinate-discretized, pseudo-spectral basis, and the Cj(t) were propagated in time 
using a fourth order variable-step Runge-Kutta algorithm, made highly efficient by 
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implementation with spatial fast fourier transform techniques. We shall use the same 
numerical methods throughout our article. 

The denumerably infinite number of conserved quantities in the ID NLSE makes 
solitons robust. Single solitons, in whatever equations they arise, are in general stable. 



Here in both the repulsive and attractive cases soliton-trains are stable |]28|. In the 
attractive case the two nodeless solution-types are cyclically stable. That is, in response 
to stochastic perturbation they evolve in a complicated manner but come asymptotically 
close to their original state in an aperiodic cycle. In figure § we illustrate the stability 
properties of the stationary states by introducing 1% stochastic white noise into all five 
fundamental stationary solution types and observing their evolution. 

Stochastic perturbation of density-notch solitons can cause a velocity shift or small- 
scale emission of phonons |4^. The phase relationships between solitons, however. 



remain intact. The depth of a density- notch is a simple function of its velocity. Therefore 
the soliton-train changes in depth but not in spacing or form. It drifts as a unit in a 
diffusive manner. This is true for both the real and complex solution-types. 

Stochastic perturbation of bright solitons can cause a velocity shift, small-scale 
emission of phonons, and/or a phase shift Although a single soliton is stable. 



multiple solitons are cyclically stable. As the phase difference between adjacent solitons 
changes from zero to vr their interaction changes from attractive to repulsive I^Sf, and it 



is this changing interaction that leads to cyclical stability. The anti-symmetric solution 
is stable. The repulsive interactions between solitons lock them into place. In the other 
two solution-types the solitons attract and repel as their relative phases are perturbed 
by the phonons that make up the white noise. We have not shown their recurrence here. 

3.2. Expectations in higher dimensionality 

We expect that in quasi- ID the five solution-types described analytically in section ^ 



and numerically in section will have the same properties as in ID. We also expect, 
and will verify in the following section, that excellent quasi-lD stationary states can be 
found as the direct product of eigenstates of the ID NLSE in each of the dimensions, 
even though the NLSE is only approximately separable. 

To reiterate our criteria for quasi- ID from section |2?1| : Ly,Lz ~ ^eff] Ly, « L^. 
In the repulsive case a non-overlapping 2-soliton solution becomes possible when the 
length is 27r^e// [10|. This provides a criterion for the breakdown of the quasi- ID 



approximation: if Ly, > Svr^e// a transverse excitation becomes likely. 

Solutions to the attractive NLSE become unstable when the negative mean field 
energy is not sufficiently opposed by the positive kinetic energy. This may be seen by 
performing a variational calculation. So one criterion we could use is the length for 
which the transverse solution changes from positive to negative energy. This is about 
2i^e//- In the literature it has been found that the 3D attractive NLSE is more likely 



to go unstable than the 2D NLSE ^J]. So we expect that the quasi-lD approximation 



will break down at a smaller transverse length in 3D than in 2D. In the attractive case 
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we expect that out of the quasi- ID regime the wavefunction should collapse. 

We note that in the repulsive case the criteria for separability and quasi-lD are not 
identical. In the limit as the transverse lengths become very large the NLSE becomes 
exactly separable for a constant transverse solution. However, the separable solutions 
can be unstable, as we shall show in section |^. 

In the repulsive case it has been shown in the mathematical literature that band 
soliton and planar soliton solutions of the NLSE in 2D and 3D can decay into vortices 
via long- wavelength transverse modulations known as snake instabilites |^5|, They 
are called band or planar solitons both to describe their shape in 2D and 3D and 
to differentiate them from ID solitons of the type formally obtained by the Inverse 
Scattering Transform [^, |3^. Thus as we increase the length of the transverse 
dimensions out of the quasi-lD regime and into 2D and 3D we expect to find that 
our soliton-train solutions break up into vortices. The BEC conserves vorticity. When 
the band or planar solitons decay the circulation quantum numbers of the resulting 
vortices should add to zero. Since the energy of a vortex is proportional to its quantum 
number squared we expect the decay to be into pairs of singly-quantized vortices of 
opposite circulation. 



4. Quasi-one-dimension 

To demonstrate the experimental viability of the ID stationary states described above 
we add one or two transverse dimensions of scale size ^e// and study the stability and 
phase engineering of solitons and stationary states in quasi-lD. The connection between 
our numerical studies and present BEC experiments is made clear. We also illustrate 
the phase rigidity of the quasi-lD BEC. We hope the following work will encourage 
experimentalists to create a quasi-lD attractive BEC and to study density-notches and 
phase rigidity in the quasi-lD regime of the repulsive BEC. 



Connection with experiments 

In table we show the length and time scales and the number of atoms for the four 
atomic species of BEC. The below conversion formula allows one to find the number of 
atoms given the healing length, trap size, and type of BEC. 

N = ^effLxLyL^ ^^^^ 
Svra 

where the healing length ^e// and the scattering length a must be in the same units and 



the trap dimensions Lx, Ly, and Lz are in units of ^e//- From section |2.1| the time scale 



is: 



2MC, 



t = ZlllAli (12) 

where M is the mass of the atomic species, t are the natural, unitless time units we use 
in the numerical studies. 
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In the table we choose ^e// = 3.0 ;um to make sohtons in the condensate 
easily observable. This is the order of the healing length used in recent soliton 



experiments |23, |T|. The time scales are on the order of 0.1 to 1 ms per natural 
time unit. Since we present results on time scales from 25 to 100 natural time units all 
our numerics take place on the experimental time scales of from 1 to 100 ms. Similiar 
results persist to 1000 natural time units. As BECs last from 1 to 75 s our studies 
extend over the lifetime of current BECs 

For computational reasons we have used somewhat smaller traps than current 
experiments, 10x1x1 or lOx^xi in the attractive case and 25x1x1 in the repulsive case. 
Thus the number of atoms is small. In the attractive case the trap can be made longer 
and remain in the quasi-lD regime as long as the number of peaks in the solutions 
is increased. In the repulsive case the trap size could be increased in each transverse 
dimension by a factor of 5 and arbitrarily in the longitudinal dimension, resulting in a 
factor of 100 or more increase in the volume. Therefore the numbers we have used in 
the simulations scale into the experimentally accessible regime. 

The assumption of separability in quasi-lD is made by projecting the total 
wavefunction onto the ground states of the transverse dimensions and integrating. This 
produces an effective healing length which depends on the transverse width: 

^ ( . 3(EK)/iyK)-i)- ^^3^ 



?>{E{mt)/K{mt)-lf 
^-ff ^-2{l + mt)E{mt)/K{mt)-{2 + mt) 

,2 _( 3{E{mt)/K{m,)-{l-mt)f 

^^ff ^ 2{-l + 2mt) E{mt) / K (mt) - {2- 5mt + 3m^t^^ ^ ^ ' 



for the repulsive and attractive cases, respectively, rut is the Jacobian elliptic parameter 
which determines the width of the transverse wavefunction, which is a sn or cn function 
described by equations (H) or (H). E{rn^) and K{mt) are complete Jacobian elliptic 
integrals. 

In the limit as -^0^, ^^jj- |.^^ in both cases and the transverse wavefunction 
is sinusoidal. As — > 1~, ^g^j — >• C,'^ and monotonically in the repulsive 

and attractive cases, respectively. For the purposes of the present quasi-lD numerical 
investigations we shall use a transverse healing length of 1 or |. For these healing lengths 
nit ~ 0.1 and ^^^^ ~ |^^. We have used these rescalings in determining the number of 
atoms in table |l|. 

Besides the ratio of the box lengths to the healing length, the numerical studies 
do not depend on the choice of the experimental parameters. In all figures we have set 
^eff = 1 and thus are working in length units scaled to the effective healing length. 

4-2. Stability 

As in section ^TT] we again add 1% stochastic white noise to the five solution-types, but 
in both the 2D and 3D quasi- ID limits. The initial state is created from the product of 
box transverse wavefunctions of the form described by equations and (^) with a box 
or periodic longitudinal wavefunction of each of the five types shown in figure |l]. 
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In figure ^ we show results of simulations in the 3D case. Trap parameters which 
enforce the quasi-lD dynamics are used: a volume of 25x1x1 in the repulsive case and 
10x|x| in the attractive case. The stability properties are identical with those found 
in ID. Although we have not shown it here, the 2D case displays essentially the same 
results. 

The repulsive solutions drift but are otherwise unchanged. The complex solution- 
type moves with a steady velocity to the left because the approximation of separability 
has failed to give the right momentum boost of the condensate to counteract the motion 
of the soliton-train. But clearly a larger momentum boost would return it to a stationary 
state. Note that the phase relationships are quite stable: constant between notches in 
figure §(a) and changing linearly between notches in figure §(c). 

The attractive anti-symmetric solution is locked into place and is stable as it was 
in ID. The two nodeless solution-types break up and recur. The symmetric, real type 
shown in figure ^(g)-(h) shows an example of this recurrence. In the third time slice 
the left-hand peak has spread out and changed its phase relationship to the right-hand 
peak. In the fourth time slice it has returned to its initial state. 

As predicted, there are small fluctuations due to the assumption of separability in 
our choice of initial state. In the repulsive case we found that beyond six healing lengths 
the solutions became unstable in both 2D and 3D. In the attractive case we found that 
the requirements for the quasi-lD regime were more complicated. A transverse length 
of \ieff was required in 3D while ^e// was sufficient in 2D. In general the stability 
depended both on the number of peaks and the solution-type. The anti-symmetric 
solution-type with a large number of peaks was the most stable configuration. This 
kept the density spread out and the contribution of the mean field term in the NLSE 
small. Other researchers have found similiar results [TO . 



Ji^.S. Manipulation via phase 



In section |2.3| we discussed simple phase profiles which impart a velocity to soliton- 
trains. Setting the whole train in motion around a thin torus creates a stable, structured 
super-current in both repulsive and attractive cases. Measuring the velocity of the 
train is a way to verify the superfiuidity of the gaseous BEG. One method would be 
combination of the new toroidal trap, created to Bose condense Cs |]23|, with optical 



phase-imprinting Note that the superfiuid properties of the attractive BEG are 

presently unknown. 

One may study soliton interactions by treating members of the train separately. 
In figure ^ we show that planar solitons in quasi- ID interact in a manner similiar to 
solitons in ID. In the repulsive case we have used an 0.37r phase jump to send two 
density notches towards one another. As may be seen in the figure, they undergo elastic 
repulsive interaction twice over 50 natural time units. 

In the attractive case we have used a real symmetric initial state with a phase ramp 
of 2ti on each bright soliton. Their first and third interactions show a characteristic 
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smearing of the density which indicates a phase difference between zero and vr p8[. In 



the second interaction they add coherently as they pass through each other, indicating 
a phase difference of zero. In their fourth interaction they repel. Therefore it is 
quite apparent that their relative phase has drifted by approximately vr, due to the 
1% stochastic noise added to the initial state. 

For computational reasons we chose phase profiles which caused interaction on a 
short time scale. But similar interactions could be observed on longer time scales by 
the appropriate choice of phase jump or ramp. Therefore solitons and stationary states 
in quasi-lD are connected in the same way as in ID. Soliton interactions can be studied 
experimentally in the BEC. 

By use of phase ramps one may also study how the BEC responds to shock. 
The repulsive case has the property of phase rigidity^ related to symmetry-breaking 
of the condensate phase. Rigidity is common to many symmetry-breaking materials 
in condensed matter physics, from liquid crystal to chalk The BEC exhibits this 
property by shattering into well-defined phase domains, the size of which are on the 
order of the healing length. In quasi- ID the borders between domains are made of 
density-notch solitons. 

To illustrate this property we use longitudinal box rather than periodic boundary 
conditions, with volume 25x1x1. A large phase ramp of = l.Gvrx is applied to a 
two- density-notch stationary state. In figure |] we show the time evolution of such a 
configuration. The first time slice shows the initial state. The second shows a shock 
wave, apparent in the density. The final time slice shows a highly excited state with 
well-defined phase domains separated by density notches. Because the kinetic energy of 
this state is very high the nonlinear term in the NLSE does not contribute signifigantly. 
Thus a linear beat phenomenon is seen in the density. We have used short time scales 
for convenience, but a smaller phase ramp would cause the same phenomenon to occur 
on a much longer, experimentally observable time scale. 



5. Beyond quasi-one-dimension 

In this section the transverse dimensions are extended beyond the quasi- ID regime in 
both two and three dimensions. Band and planar solitons become unstable, and shocking 
the condensate produces chaos. 



5.1. Vortex creation 

In figures and ^ we show the decay of band and planar solitons into vortices. In both 
cases we begin with 2-soliton-train stationary states made in the same way as described 
in section ^ We use dimensions 25x12 in 2D and 25x12x12 in 3D with periodic boundary 
conditions in the longitudinal direction. White noise causes the solitons to decay into 
a pair of vortices of opposite circulation, conserving vorticity. The vortices are singly 
quantized, as may be seen in the figures by traversing a node in the phase. The colors 
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pass once around the phase color circle, which is defined to be 27r in extent. In the 
plots we chose times that would highlight the decay. However, the vortices are stable, 
spatially stationary, and do not interact over the full length of the study. 

White noise in the BEC is principally due to uncondensed atoms and the number 
of such atoms is a function of temperature. Thus noise is nominally an experimental 
parameter. The time from creation of the initial stationary solitons to their decay into 
vortices is a function of the noise. With the correct choice of temperature and volume. 



stationary band or planar solitons, made via phase engineering can be used to 
create vortices in a controlled manner. 



5.2. Vortex interaction 

Band and planar solitons respond to simple phase profiles in the same way as solitons 
in quasi- ID. The following method permits direct study of vortex interactions. 

As explained above, the time from creation of the initial state to decay into vortices 
is a function of the temperature of the BEC. As in quasi-lD, they move toward each 
other with a velocity which depends on the phase jump. Thus two solitons can be made 
to decay into vortices shortly before collision. An example of this process in 2D is shown 
in figure ^ We use a phase jump of 7r/15, a noise level of 0.1%, and a container identical 
to the one in figure |^. In the first time slice a snake instability and a vortex pair have 
developed. By t = 36 the vortices are strongly interacting. They have become a pair of 
vortex dipoles. By t = 48 they have completed their interaction and are moving away 
from each other around the torus. The time slices in the figure are evenly spaced. Thus 
it can be seen that, unlike in quasi-lD, the collision is not elastic. The relative vortex 
velocity in the first three time slices, before the collision, is not the same as in the last 
two. In fact the two vortex dipoles no longer have the same speed. 

Although we have not pictured it here, the evolution of the above collision is very 
intriguing. Despite the complicated density profile, it is apparent from the phase that 
the remnants of the two band solitons oscillate between single vortices, vortex dipoles, 
and band solitons. The relationship between the spectra of these various collective 
excitations is an unsolved problem. 



5. 3. Shock and phase rigidity 

Explorations of 2D and 3D out of the quasi- ID regime have thus far have been with 
periodic boundary conditions in the longitudinal direction. If instead box boundary 
conditions are used, together with a phase ramp of the kind desribed in the previous 
section, the BEC becomes chaotic. As in quasi-lD, it exhibits the property of phase 
rigidity. In response to shock it breaks up into well-defined phase domains. 

To illustrate this effect a volume of 25x12 with box boundary conditions in both 
dimensions was used. A large phase ramp of 4>{x,y) = l.Gnx + 0.27iy pushed a two- 
soliton, approximate stationary state into the walls at a skewed angle. The time 



evolution of the resulting state is shown in figure 10. The first time slice shows the 
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initial state. The second time slice shows a shock wave. By the final two time slices the 
condensate appears to have fractured into a chaotic assortment of phase domains. The 
phase appears to be completely decoherent at scales larger than the healing length. 



6. Conclusion 



We have analytically described and numerically illustrated the experimental viability 
of soliton-train stationary states in quasi-one-dimension. Soliton-trains can be 
manipulated with simple phase profiles. In this way one can create a structured super- 
current or study soliton interactions. 

In the repulsive BEC quasi- ID soliton-train solutions are stable. In the attractive 
BEC those in which the members of the train are separated by nodes are stable. Other 
solution-types are cyclically stable. The attractive BEC is stable against collapse in the 
quasi- ID regime. 

Beyond the quasi- ID regime we have shown that band and planar solitons decay 
into vortices when perturbed by white noise. This is not only a new method of vortex 
creation but also allows for the direct study of vortex interactions. We have also shown 
that the phase of the BEC is rigid; when shocked it shatters chaotically. 
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Figure 1. Shown are the five stationary solution- types in one dimension under periodic 
boundary conditions. Repulsive case: (a) amplitude of real sohition and (b) amplitude 
and phase of complex solution. Attractive case: (c) amplitude of real symmetric and 
antisymmetric solutions and (d) amplitude and phase of complex solution. All plots 
are normalized to the unit interval. 
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(a) (b) 
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Figure 2. The five stationary solution- types in one dimension under periodic 
boundary conditions are propagated in time numerically for 100 natural time units. 
Shown are phase and density for (a)-(b) real solutions, repulsive case; (c)-(d) complex 
solutions, repulsive case; (e)-(f) anti-symmetric real solutions, attractive case; (g)-(h) 
symmetric real solutions, attractive case; (i)-(j) complex solutions, attractive case. 





Figure 3. The five stationary solution- types in one dimension under periodic 
boundary conditions arc propagated in time numerically, as in figure ^, for 100 natural 
time units. Stochastic noise at the level of 1% is added to the initial state, (a)- 
(f) are stable; (g)-(j) are cychcally stable. Shown are phase and density for (a)-(b) 
real solutions, repulsive case; (c)-(d) complex solutions, repulsive case; (e)-(f) anti- 
symmetric real solutions, attractive case; (g)-(h) symmetric real solutions, attractive 
case; (i)-(j) complex solutions, attractive case. 
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Figure 4. The five stationary solution-types are extended into quasi-one-dimension 
under periodic boundary conditions. They are propagated in time numerically as 
in figure |[ for 50 and 25 natural time units in the repulsive and attractive cases. 
Stochastic noise at the level of 1% is added to the initial state, (a)-(f) are stable; (g)- 
(j) are cyclically stable. Shown are four time slices through the xy plane of phase and 
density of fully 3D simulations for (a)-(b) real solutions, repulsive case; (c)-(d) complex 
solutions, repulsive case; (e)-(f) anti-symmetric real solutions, attractive case; (g)-(h) 
symmetric real solutions, attractive case; (i)-(j) complex solutions, attractive case. 
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Figure 5. In quasi-lD soliton-train solutions can be manipulated by simple phase 
profiles. Here we show examples of such manipulation for (a)-(b) the repulsive and 
(c)-(d) the attractive cases. The solitons have been given equal and opposite velocity. 
The evolution of their phase and density is shown by plotting the xy plane at the 
midpoint of the z co-ordinate in time slices running down the page. We have added 
1% stochastic noise to the initial state. The results shown follow from a fully 3D 
propagation with dimensions 25x1x1 in the repulsive case and lOx^x^ in the attractive 
case. 
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Figure 6. A rapidly translating condensate in a quasi-lD box is suddenly stopped. 
The condensate shatters into domains of constant phase due to the fragility of this 
phase-symmetry-breaking system. 
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Figure 7. In 2D band solitons decay into pairs of oppositely charged vortices. We 
have added 0.01% stochastic noise to an initial two-soliton stationary solution in order 
to demonstrate this effect. Shown are five slices of the xy plane at the midpoint of the z 
co-ordinate, equally spaced in time. The longitudinal direction has periodic boundary 
conditions. Note that the phase leads the instability. 




Figure 8. In 3D planar solitons decay via a snake instability. We have used the same 
noise level as in figure |^. We have highlighted an earlier set of times, as the decay 
occurs earlier than in 2D. Shown is the xy plane at the midpoint of the z co-ordinate. 




Figure 9. Collision of two initial band solitons which have evolved into two oppositely 
charged vortex pairs. We have added 0.1% stochastic noise. As in figure ||the relative 
velocity is induced by phase discontinuities across the band solitons, but the interaction 
is that of two vortex dipoles. This is an example of how to use the noise-induced 
instability of band solitons in 2 and 3D to study vortex interactions. 
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Figure 10. A rapidly translating quasi-2D repulsive condensate is stopped. Following 
propagation of a density-dependent shock front, a fully chaotic phase and spatial 
distribution is obtained. 



